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^ '. Abstract 

I We describe a scheme into which a camera is turned into an efficient tunable frequency filter of a 

' jy^" few Hertz bandwidth in an off-axis, heterodyne optical mixing configuration, enabling to perform 

o ■ 

parallel, high-resolution coherent spectral imaging. This approach is made possible through the 

^: 

O , combination of a spatial and temporal modulation of the signal to reject noise contributions. Ex- 

O ' perimental data obtained with dynamically scattered light by a suspension of particles in brownian 

J>-j' motion is interpreted. 
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INTRODUCTION 



Coherent spectral imaging 



Coherent spectroscopy enables one to study mechanisms involving dynamic light 
scattering. The spectral distribution of a monochromatic optical field scattered by moving 
particles is modified as a consequence of momentum transfer (Doppler broadening). The 
measurement of the Doppler linewidth of this field (referred as object field) with an optimal 
sensitivity is crucial, since Doppler conversion yields are typically low. 



Optical mixing (or postdetection filtering) techniques are derived from RF spectroscopy 
techniques They can be grouped in two categories 0, 3] : homodyne and heterodyne 
schemes. In homodyne mixing (fig. [U^a)), self-beating object light impinges on a Aujpu- 
bandwidth photodetector (also referred as optical mixer or photo- mixer, PM). To assess 
a frequency component of the object light, the output of the PM is sent to a spectrum 
analyser, whose bandwidth Au-p defines the detection resolution. The resulting spectra 
are proportional to the second-order object field spectral distribution S2{uj). In 

heterodyne mixing, sketched in fig. [I](b), the object light is mixed onto a PM with a 
frequency-shifted reference beam, also called local oscillator (LO). The LO field (-Elo) is 
detuned to provoke a heterodyne beat of the object-LO field cross contributions to the 
recorded intensity. This beat is sampled by a PM and sent to a spectrum analyser whose 
bandwidth Aujp defines the apparatus resolution. This scheme enables one to measure the 
first order object field spectral distribution Si{uj) 



In heterodyne optical mixing experiments, the PM bandwidth AtupM defines the span of 
the measurable spectrum (usually ~ 1 GHz). The resolution Aup can be lowered down to 
the sub Hertz range, which is suitable for, amongother applications, spectroscopy of liquid 
and solid surfaces [41], dynamic light scattering 21 and in vivo laser Doppler anemometry 
jsl. But these schemes are inadequate for imaging applications, because measurements are 
done on one point. The absence of spatial resolution has been sidestepped by scanning 
techniques 0, Q], designed at the expense of temporal resolution. 



2 



We present a heterodyne optical mixing detection scheme on an array detector (configu- 
ration of fig. [11(b)). Typical array detectors sampling rates seldom run over 1 kHz, failing 
to provide a bandwidth large enough for most Doppler applications to date. Nevertheless, 
their strong advantage is to perform a parallel detection over a large number of pixels. We 
present a spatial and temporal modulation scheme (spatiotemporal heterodyning) that uses 
the spatial sampling capabilities of an area detector to counterbalance the noise issue of a 
measurement in heterodyne configuration performed in the low temporal frequency range 
(e.g. with a typical CCD camera). The issue of using narrow-bandwidth camera PMs is 
alleviated by detuning the LO field optical frequency accordingly to the desired spectral 
point of the object field to measure. Post-detection filtering results from a numerical 
Fourier transform over a limited number of acquired images. 



This heterodyne optical mixing scheme on a low frame rate array detector can be used 
as a filter to analyze coherent light. It has already been used in several applications yet. 



ing_th( 

isM, 



includin g th e detection of ultrasound- modulated diffuse photons 
analysis 



low-light spectrum 



study 



laser Doppler imaging [ll, 1^, isj, and dynamic coherent backscattering effect 



14| . The purpose of this paper is to present its mechanism. 



Time domain description of the fields 

We consider the spatially and temporally coherent light field of a CW, single axial mode 
laser (dimensionless scalar representation) : 

E^it) = ALA(t) exp {iuLt + i0(t)) (1) 

where cul is the angular optical frequency, is the amplitude of the field (positive constant), 
At(t) = 1 -|- a(t); \a(t)\ ^ 1 describes the laser amplitude fiuctuations and 0(t) the phase 
fiuctuations. This field shines a collection of scatterers that re-emits the object field (or 
scattered field), described by the following function : 

Eoit) = AoXt(t) A(t) exp {iuji^t + i^it)) (2) 

where Aq is a positive constant and Xt{t) is the time-domain phase and amplitude fiuc- 
tuation induced by dynamic scattering of the laser field, i.e. the cause of the object field 
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fluctuations due to dynamic scattering we intend to study. As in conventional heterodyne 
detection schemes, a part of the laser field, taken-out from the main beam constitutes the 
reference (LO) field : 

ELo{t) = Aj^oAtit) exp (zc^lo^ + i<P{t)) (3) 

where Alo is a positive constant. The LO optical frequency is shifted with respect to the 
main laser beam by culo — to provoke a tunable temporal modulation of the interference 
pattern resulting from the mix of the scattered and reference fields. It can be done 
experimentally by diffracting the reference beam with RF-driven Bragg cells (acousto-optic 
modulators) for example. 

The expression of the light instant intensity impinging on the camera PMs is : 

t{t) ^ toiRe[E])^ (4) 

where io is a positive constant and Re[i?] is the real part of the relevant optical complex field 
E. The PMs are considered point-like, their antenna properties are out of the scope of 
this paper. To take into account the frequency filtering implication of finite time-domain 
integration, the average intensity detected by the square-law camera PMs is : 

I{t) = lo {E{t + T) + E*{t + r))2 dr (5) 

where Iq is a positive constant and Tg the exposure time. In the frequency- domain, this 
integration corresponds to a sinc-shaped low-pass filter function whose bandwidth (defined 
as the distance between the peak and the first zero) is l/rg. 



OPTICAL CONFIGURATION 



Setup 

The common name for optical mixing experiments on an array detector is digital 
holography. The original underlying interference technique was invented to improve 
electron microscopy resolution and was successfully applied to monochromatic light imaging 
[l^. Since the availability of lasers and then digital cameras, digital holography [l7| has 
become an integrant part of many coherent imaging schemes, for its propensity to record 
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both quadratures of an optical field. 



The optical confieruration we use can be either a lensless Fourier off-axis holography setup 

nn ^ 

[18|, Il9| or a Fresnel off-axis holography setup [20]. In the lensless setup, sketched in fig. |21 

and considered throughout this study, the point source of the spherical reference wave (LO 

focal point) is located at (xg,?/g) in the plane of the object to match the average curvature 

of the object field. An interference pattern is recorded in the camera plane {x,y). The 

reconstruction algorithm used to calculate the spatial distribution of the scattered field in 

the object plane consists of only one spatial fast Fourier transform (FFT) [3]. This setup 

is used to avoid the need of a general Fresnel reconstruction algorithm. In the Fresnel setup 
□ 

[20|, the LO is a plane wave and the image reconstruction requires the use of numerical 
lenses and at least two successive spatial domain FFTs 21|. 



The digital hologram is recorded with the frequency-shifting method introduced in the 
heterodyne holography technique 22], used in several imaging and detection schemes 
10, [3] for its ability to discriminate efficiently a few Doppler-shifted photons, according to 
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their frequency, from background light [8|, llO| . The frequency-shifting method was chosen 
as an accurate dynamic phase-shifting detection scheme. This method is derived from 
static phase-shifting interferometry [23] applied to phase- shifting holography [24]. It consists 
of detuning the LO optical frequency with respect to the main laser beam, to provoke a time 
domain modulation of the intensity pattern resulting from the cross contribution of object 
and LO fields. 



Spatial and temporal dependencies 

Space (subscript s) and time (subscript t) dependencies of X and A are noted X{x, y, t) = 
Xs{x,y) ■ Xt(t), A{x,y,t) = As{x,y) ■ At{t). Spatial frequencies (subscript sf) and temporal 
frequencies (subscript f) reciprocal distributions form Fourier pairs with time domain and 
spatial domain distributions. We have : 

l-T/2 

Y{{uj)= Yt{t) exp{-iujt) dt (6) 

and 
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X exp(— 227r(x^ + yr])) dxdy (7) 

where Y is the considered distribution, T is the measurement time, Ax and Ay are the 
camera sensor widths. 

In the lensless Fourier configuration, the LO wave curvature matches the average curva- 
ture of the object field. Therefore, the heterodyne interference pattern on the detector is 
formally equivalent [18] to the one that would impinge on the detector under Fraunhofer con- 
ditions. Consequently, we can consider the (far field diffraction) equivalent situation where 
the LO wave is a tilted plane wave in the detector plane and the object field distribution 
is the Fourier transform of its spatial distribution in the object plane. The object and LO 
complex optical fields take the following form: 

Eo{x,y,t) = 
AoXs{x,y)As{x,y) 
xXt(t)At(t)exp(zcuLt + «0(t)) (8) 

ELo{x,y,t) = 
Ai^oAix, y) exp (2m(x^o + VVo)) 

xA(t)exp(itULot + (9) 

where = ^o/i^d) and r/o = Vo/i^^)- The tilt angle of the reference field leads to the phase 
factor exp (2i7r(x^o + Wo))- The detected intensity is : 

I{x,y,t) = 

h / [^0(2;, i + ^) + Ei^o{x, y,t + T) 

+E*o{x, y,t + T) + Eloix, y, t + r)]^ dr (10) 

Expression [10] has 16 terms among which 10 vanish because of the presence of optical fre- 
quencies in the phase factors j3] , which lead to oscillations outside the detection bandwidth. 
It can be rewritten as : 



I{x,y,t)/lQ = 
AoAi^oh{x,y) ■ h{t) 
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+AoAi^ohix,y) ■ hit) 
+AloIsix,y)-hit) 

+Alh{x,y)-h{t) (11) 



where the time-domain fluctuations of the intensity are : 

rTo/2 



hit) = r xt{t+T)\At{t+T)f 

J~Tc/2 

X exp (-z(a;LO - t^L)(i + r)) dr (12) 
m= r^'x:{t + T)\A{t + T)f 

"'-re/2 

X exp (i(u;Lo - ^^L)(t + r)) dr (13) 



hit)= |A(t + r)|'dr (14) 



h{t)= \X,{t + r)f\A,{t + r)fdT (15) 

J-Tc/2 



rrc/2 

and the spatial contributions to the intensity take the following form 



h{x,y) = 

I A(x, y) f Xs(x, y) exp (2z7r(a;^o + Wo)) (16) 

h{x,y) = 

\A{x,y)f X*{x,y) exp {-2i-iT{x^o + yr]o)) (17) 
where = Xq/^XcI) and ?7o = yy{\d). 

h{x,y) = \A,{x,y)f (18) 

h{x,y) = \X,{x,y)\^\A,{x,y)\^ (19) 

Temporal and spatial dependencies can be treated separately. Eq. [TT] has four terms : the 
first two are object-LO field cross terms. The third and fourth terms correspond to the LO 
and object fields self beating contributions. The two first spatial contributions, h{x,y) & 
h{x, y) will be referred as heterodyne, and h{x, y) & h{x, y) as homodyne, as well as their 



temporal counterparts. 

n 

An assumption about the relative amplitudes of the fields is usually made |2i] : the object 
field amplitude is much lower than the LO field amplitude (heterodyne regime) : 

< ^Lo (20) 

This assumption leads to neglecting the object field self-beating contribution, because its 
relative amplitude with respect to the cross terms is Aq/Alo ^ 1- Additionally, in usual 
optical mixing schemes on photodiodes, laser amplitude fiuctuations are neglected (i.e. 
a{t) = 0) Q], since the measurement is done with detectors whose bandwidth is large 
enough to make this assumption valid, as the detection is made at higher frequencies. But 
with a slow PM such as a CCD camera, the frequency- domain noise resulting from these 
fiuctuations is not negligible. 

A sequence of n images is sampled and leads to a data cube I{x, y, t), acquired for a given 
detuning frequency culo — ^l- To make a spatial map of one frequency component of the 
object field, a demodulation in the time and spatial domains is performed. As we will see, 
the heterodyne holography scheme allows encoding of the spectral and spatial information 
about the object field in the reciprocal space of the data cube, and it also allows spatial 
discrimination of the self beating terms from the object-LO field cross terms. 

SIGNAL DEMODULATION IN THE TIME DOMAIN. 

Ji and I2 are referred as the heterodyne terms (modulated at the detuning frequency 
cjlo ~ ^l, on purpose). They correspond to the object and LO field cross terms. I3 
and J4 are referred as the homodyne terms (not explicitly modulated at the detuning 
frequency). These terms correspond to the object field (and the LO field) self beating 
intensity contributions. 

To measure a spectral component of the object field, n samples of I{x,y,t) along the 
time axis are acquired. These images are used to calculate both quadratures of the object 
field by making a time domain demodulation consisting of calculating the first harmonic 
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component {I {x , y , us/ n)) of the recorded sequence. This method is also used in temporal 
heterodyne inline holography 25|, l26|] . The relation between the result of this demodulation 
and the temporal frequency content of the object field is investigated. 



First heterodyne term 

Let's consider the first heterodyne term : 

h{t) = F{t) exp(-zKo - ^l)^) (21) 

where F{t) is : 

F{t)= \A,{t + T)\^X,{t + T) 

X exp(-i(u;Lo - dr (22) 

According to eq. [221 F{t) corresponds to the result of the application of a bandpass filter 
(of center frequency culo — and bandwidth 1/re) on Xt{t). 

To demodulate the signal, the first harmonic {uJs/n) of the sampled image sequence is 
calculated by FFT : 

n 

hiujs/n) = J2 Htk) exp{-2tkTT/n) (23) 
fc=i 

where = 2fc7r/c<Js is the instant at which the k^^ image is recorded. We have : 

n 

hi^^s/n) = F{tk) exp(-i(t^LO - tUL)tfc) 

k=l 

X exp(— 2iA;7r/n) (24) 

Ii{uJs/n) is a measurement of the culo — + ^^s/n frequency component of F(t) (the corre- 
sponding instrumental bandwidth is uj^/n). Eq. [221 and [2H show that the spectral resolution 
of this measurement should either be limited by the measurement time (n/us) and/or the 
spectral distribution of the amplitude noise of the laser. 



Second heterodyne term 

The assessment of the quantity l2{uJs/n) is straightforward : 

n 

hiujs/n) = F*(tk) exp(2(u;LO - cul)^) 

k=l 
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X exp(— 2ifc7r/n) 



(25) 



where F*{t) is the complex conjugate of F{t) defined by eq. [221 Equations [221 and [251 define 
the action of two successive selective filters on Xt{t). The filters bandwidths are the same 
as the ones introduced to describe the first heterodyne term. Their action is presented and 
compared to their counterparts appearing in the first homodyne term in the next section. 

Bandpass filtering 

If we neglect the contribution of the amplitude fiuctuations At{t) to the laser linewidth, 
the impulse response functions for the first and second heterodyne terms {B^{uJi^o — cul) 
and B^{ui^o — c<jl) respectively) calculated by setting Xt{t) = 1 in expressions [2^ [251 and 
[22ltake the form : 

B±{uJi,o - ^^l) = sine ((culo - c^l) ^o) 
X ^ exp(— 2iA;7r/n) exp(=F2iA;7r^^^^^^^^^ — —) (26) 

k=l 

The instrumental response \B±{uJi^o ^ ^hjf was measured experimentally, by analyzing 
light backscattered by a static object. Data points of this instrumental response are 
plotted as a function of the detuning frequency in fig. [3l The square amplitude of the 
function described by eq. [261 is represented in fig. [H Although their shape is similar, the 
dynamic range of the theoretical filter is wider than the measured response. The difference 
is attributed to the linewidth contribution of laser intensity fluctuations. One important 
thing to remark about those instrumental responses is their dissymmetry. 

\Ii{us/n)\^ and |/2(c<Js/n)p are calculated to represent quantities homogenous to the 
optical power of the object fleld. This power is proportional to the object fleld power 
resulting from the integration of its spectral density in the B± windows. We call these 
distributions signal (or true image) and ghost (or twin image), respectively. B± describe 
bandpass fllters of width Aup = cus/n, centered on culo — ^^l ± ^^s Z''^- Consequently, this 
scheme allows one to measure the t^Lo — i^l ± t^s/""- frequency components of the object 
fleld fluctuations Xt{t). 
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If we set the detuning frequency to : 

t^LO - = Au; - us/n (27) 

the signal will correspond to the Auo frequency component of the field and the ghost area 
to the Atu + 2uJs/n component. We have \Ii{ijjQ,/n)\'^ ^ 5*1 (Acj) and |/2(co's/n)p ^ S'i(Aa; + 
2uj^/n). The corresponding frequency diagram is sketched in fig. [51 

Homodyne terms. LO and object fields self beating 

In usual heterodyne monodetection schemes, self-beating contributions are considered as 
an unavoidable noise component, which can actually be neglected in a high bandwidth PM 
detection scheme. When using low frequency PMs such as camera pixels, these homodyne 
contributions (defined by eqs. [H] and [T5l) have to be taken into account. We show in the 
next section how these terms are filtered-out spatially in the detection process. 

SIGNAL DEMODULATION IN THE SPATIAL DOMAIN. FILTERING OF SELF- 
BEATING CONTRIBUTIONS 

Heterodyne contributions 

The precious advantage of holography is to allow sampling of the diffracted object 
complex field (in phase and amplitude), whose spatial distribution in the detector plane is 
described by the X^^x, y) function. The time domain demodulation presented in the previous 
section enables one to measure the spatial distribution of a tunable frequency component 
of the object field, in amplitude and phase (i.e. both quadratures). Thus, after tem- 
poral demodulation, the actual complex distributions defined by eqs. [12] and [T7] are available. 

As a result of the off-axis configuration, the object-LO field cross terms (or heterodyne 
contributions) carry phase factors exp (±2i7r(a;^o + Z/^o))- The relative distance between 
each point {x',y') of the object (described by Xs{x',y')) and the LO focal point (xQ,yQ) in 
the object plane is encoded in a set of parallel (complex) Young's fringes in the detector 
plane (x, y). The periods of these fringes are (Ad/ (x' — x'q) and \d/{y' — Vq))- Their spatial 
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frequencies are noted A^o = {x' — XQ)/{\d) and A?7o = {y' — ?/q)/(Ac?). 

In Fraunhofer conditions, the relative distances between each point of the object {x', y') 
and the LO focal point (xq,?/q) are proportional to the spatial frequencies of the hologram 
in the camera plane {^,ri). We have : 

|X,f(e = x'/iXd),r] = y'/iXd))\' cc \X,ix',y')\' (28) 

ience the requirement for only one spatial Fourier transform to reconstruct the image 

3, 
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27l | , i.e. to calculate the spatial distribution \Xs{x',y')\ . 



The laser amplitude noise leads to flat-field fluctuations in the camera plane. The spatial 
frequencies distribution of this noise is considered to be a centered dirac : 

rAx/2 i-Ay/2 ^ 

\Aix,y)\ 



-Ax/2 J~Ay/2 

X exp (^'~^)^(^-^' + yy')^ dxdy 

^6{x',y') (29) 

this assumption is implicit in lensless digital holography. The spatial frequency content of 
Ii{x,y) , assessed by a FFT, takes the following form : 

\hi^ = x'/iXd),v = y'/m)\'^ 

|Xs(x'-Xo,y'-yo)l^ (30) 

The ghost distribution is : 

\i,{^ = x'/iXd),v = y'/m)\'^ 

|Xs(xo -x',?/o -y')l^ (31) 

In the lensless configuration, both signal and ghost spatial distributions are focused in the 
same reconstruction plane. Eqs. [30] and EH define two spatial distributions flipped one with 
respect to the other and shifted away by ^{x'q, y'^) from the center of the reconstructed 
hologram. The spatial instrumental response width corresponds to less than the distance 



between two adjacent pixels in the reconstructed image 



27|; its contribution is neglected. 
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Homodyne contributions 



The object and LO self beating contributions I3 and I4 are not encoded in the spatial 
complex fringe system which results from the interference of the object field and the off- 
axis LO. Under the assumption of a perfect fiat-field amplitude noise (eq. [22]), the spatial 
frequency content of eqs. [18] and [19] take the form : 

= x'/{Xd),v = y'/m) oc 6ix', y') (32) 

and 

hi^ = x'/iXd),r] = y'/iXd))^ 

X,{x',y')*X:ix',y') (33) 

These terms are restituted inline : they are centered on {x' = 0,y' = 0) in the hologram 
reconstructed in the object plane. This enables one to discriminate them spatially from 
off-axis heterodyne contributions. 



EXPERIMENTS 



iil,Q,[i? 



In previous publications [11, ll2,llJ, the spatiotemporal heterodyne detection technique 
was used to perform parallel imaging. Here, we will focus our attention on the temporal 
domain. We will study in particular how measured frequency spectra are affected by the 
camera frame rate, exposure time, and total measurement time. 



Setup and data acquisition details 
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he experimental setup, which is sketched in FigjS] has been described previously 10|, 
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13l |. The light source is a Sanyo DL-7140-201 diode laser (A = 780 nm, 50 mW 
for 95 mA of current), the camera is a PCO Pixelfiy digital CCD camera (12 bit, frame 
rate u;s/(27r) ~ 12.5 Hz, exposure time Te < 2'k/ujs = 80 ms, with 1280 x 1024 pixels of 
6.7 X 6.7 /im), and the frequency shifter is a set of two acousto-optic modulators AOMl and 
A0M2 (Crystal Technology; ti;AOMi,2 — 80 MHz). A neutral density is used to control the 
LO beam intensity. The sample is a 1 cm x 1 cm x 5 cm rectangular PMMA transparent 
cell filled with a diluted suspension of latex spheres in water (Polybead: Polyscience Inc., 
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diameter 0.48 /im, undiluted concentration: 2.62% solids-latex). The Polybead suspension 
is diluted by a factor ^ 9 (0.5 ml of undiluted suspension + 4 ml of water). A rectangular 
aperture (7 mm x 3 cm) located just in front of the cell removes the parasitic light diffused 
along cell sides. The aperture delimiting the imaged side of the cell is located at a distance 
d ~ 39 cm of the camera. This sample is observed in transmission and the spectrum of 
the light dynamically scattered by the latex spheres in brownian motion is measured by 
sweeping the AOMl frequency so that the detuning frequency cjlo ~ "^l is swept from 
to 16 kHz in 81 frequency points (200 Hz step). For each frequency point a sequence of 
32 consecutive CCD frames is recorded to the PC computer hard disk. For each frequency 
sweep, 32 X 81 = 2592 images are recorded. 

Data analysis 

A first experiment consists of acquiring data with 3 frequency sweeps with Te = 80 ms, 
Te = 20 ms and Xe = 5 ms. For each frequency point, a n = 4 phase demodulation is 
performed by using either 4, 8, 16 or the complete set of 32 recorded images. When more 
than 4 images are used, a first n = 4 phase demodulation is done with images 1 to 4, which 
is followed by a second demodulation with images 5 to 8, and so on. The complex signal 
resulting from the demodulation of each set of 4 images is then summed to get the final 
demodulation result. The total measurement time of one frequency point vary from 320 ms 
(4 X 80 ms) to 2.6 s (32 x 80 ms), while the total exposure time vary from 20 ms (4 x Tg 
with Te = 5 ms) to 2.6 s (32 x with = 80 ms). 

Fig. [6] represents the resulting images of the rectangular aperture, illuminated in 
transmission via the diffusing cell, for an exposure time t^, = 80 ms, and for a detuning 
frequency (colo - uJiM{2t:) of Hz (a), 400 Hz (b), 4000 Hz (c) and 8000 Hz (d). The 
images are k-space images obtained by FFT of the hologram measured in the CCD plane 
and n = 4 phase demodulation with 4 images per frequency point. They are displayed 
in arbitrary logarithmic scale. The true image (signal) of the rectangular aperture is the 
region 1 in fig. [6]^a). This distribution is the spectral component of the object field of 
frequency cul ~ ^s/^- The twin image (ghost) of the aperture is the region 3. In accordance 
with eq. [30] and [3ll it is symmetric with respect to the center (tag 2) of the k-space plane 
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(null spatial frequency). The twin image corresponds to the uji^ + ujs/n spectral component. 
Increasing the detuning frequency cjlo — ^l, the true and twin images of the aperture 
become darker and darker. Nevertheless, they are still visible for a 8000 Hz offset. One can 
notice that for null detuning (fig. (Hl^a)) some statically scattered parasitic light is detected 
out of the aperture image (region 4). This parasitic light is no more detected when the 
frequency offset is non zero (images (b) to (d)). 

To perform a quantitative analysis of the signal, we have summed the signal intensity 
along the columns of the images represented on Figl6], from line 312 to line 712 (horizontal 
dashed lines on FigE](a)). These traces are represented on fig. [71 The true image corresponds 
to region 1, the twin image to region 2. Regions 3 and 4 correspond to the background 
signal. As mentioned in [a this background signal is due to the shot-noise on the LO. 
Taking into account the heterodyne gain, the background signal corresponds to an optical 
signal of one photo electron per pixel for the whole measurement sequence. The background 
signal provides here a very simple absolute calibration for the optical signal diffused by the 
sample. Here, in the center of the aperture, the sample diffused signal is about 3.5 times 
the background i.e. 3.5 photo electron per pixel. 

Spectra of Doppler-broadened light in the kHz range 

From the reconstructed images (Fig 16]), we have calculated the shape of the first order 
spectrum of the diffused object field. To get these spectra, we have calculated the sum of 
the area of the true image and twin image areas (regions 1 and 2 in the traces of fig. 17D, 
and subtracted the background (regions 3 and 4 in fig. 17]). Fig. 18] and Fig. 19] show brownian 
frequency spectra of the light dynamically scattered through the cell. These curves have 
been normalized by the area under the raw lineshapes. 

Data was collected for detuning frequencies up to 16 kHz, much larger than the hetero- 
dyne receiver bandwidth (0.4 Hz for curve (d), reciprocal of the 2.6 s measurement time). 
The shape of the spectrum does not depend on the number of images used to perform the 
demodulation, as we can see in Fig. 18] In the range of exposure times (of one frame) used 
for the measurement, the shape of the spectra reported on fig. l9]does vary with for the 
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following values of the latter : 5 ms (curve c), 20 ms (curve b), and 80 ms (curve a). FiglS] 
and Figjni which exhibit the same lineshape, show that these spectra do not depend on the 
total exposure time, which varies from 20 ms (Fig. M^c)) to 2.6 s (Fig. [H](d)). One must 
notice that these results are valid because the width of the spectrum (~ 2.5 kHz half width) 
is large compared to the instrumental response. In the all the presented results, the in- 
strumental response is always smaller than l/xe < 200 Hz (width of the sine factor in Eq j26!l . 



Additionally, we measured the first order spectrum of the object field in fixed instrumental 
conditions (4-phase demodulation, exposure time Te = 80 ms) to study different samples. 
The scattering parameters of the suspension were changed in the following manner : the 
latex beads concentration was halved from one measurement to another. The four Doppler 
lineshapes reported on fig. [10] correspond to latex volumic fractions ranging from 2.9 x 10~^ 
(curve a) to 3.6 x 10^^ (curve d). 



DISCUSSION 



The presented technique uses and goes beyond the concept of time-average holography. 
In time-averaged holography, the detector, whose exposure time is large, behaves as a 
low-pass filter of bandwidth I/tc. It selects the field components Eo{u) whose frequencies 
uj = uji^ + Auj are close to c^lO) i-e. which satisfy lAcuIxe < 1). This low-pass filtering effect is 
used to study vibrating objects like musical instruments j28| by imaging the regions which 
are not moving or which move with a given phase and velocity. In the first case (selection 
of the non moving regions), the LO beam is not shifted in frequency 29], [s^. In the second 
case (moving regions), the LO beam is modulated at the object vibration frequency 3ll. |32| 
(which, in the temporal frequency domain, corresponds to the generation of adequate LO 
frequency sidebands. 



Our scheme enables detection of a tunable frequency component of the light with a sharp 
bandwidth defined by the inverse of the acquisition time of a sequence of n images. But 
reading-out an optical beat at typically low camera frame rates is difficult in practice in 
the case of a weak-amplitude object field, because the low frequency part of the temporal 
frequency spectrum is highly noisy : it contains LO and object field self-beating contribu- 
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tions. To discriminate the signal from these noise contributions, a lensless Fourier off-axis 
configuration is used to encode, in a Young fringe system, the distribution of relative posi- 
tions of object points with respect to the reference point. This scheme allows to restitute 
the object field spatial distribution in the Fourier reciprocal space of the detection plane. 
Having recourse to an off-axis configuration to sample the heterodyne optical beat enables 
spatial frequency discrimination of the heterodyne terms (object and reference light cross 
terms) carrying useful information from the object and reference light self-beating terms 
carrying what turns out to be unwanted information. Time domain and spatial domain 
modulation/demodulation schemes present strong similarities. Equation [16] is the spatial 
counterpart of eq. Signal and ghost distributions are beating at icus/'"' in time domain. 
This beating enables to perform a ra-phase temporal demodulation of the recorded image 
sequence to assess object and ghost fields in quadrature. Since these complex fields are 
spatially modulated because of the off-axis interferometry configuration (according to eqs. 
[16] and [17]) , their distributions in the object plane are translated by ±(xo,?/o) (respectively) 
away from (0, 0)-centered homodyne contributions. This is highly valuable since it allows 
one to filter them out spatially. 



CONCLUSION 



The presented scheme, based on heterodyne optical mixing onto a parallel detector used 
as a selective frequency filter, is particulary suited to high resolution spectral imaging. 
An tunable frequency component of the object field is acquired at a time. The camera 
finite exposure time and the finite measurement time lead to a frequency domain bandpass 
filter, which width is the reciprocal of the measurement time. The small dissymmetry of 
the frequency domain instrumental response should be taken into account in quantitative 
measurements. The available range of frequency shifts at which measurements can be 
made is not limited to the photo-mixer bandwidth (contrary to wide-field laser Doppler 



331), and benefits from heterodyne amplification. The measurement is sensitive : thanks 



to the off-axis holographic setup, the self beating intensity contributions resulting from 
laser instabilities and scattered light self interference within the camera bandwidth can 
be efficiently filtered-out. The propensity to separate spatially cross terms intensity 
contributions from self beating contributions lies in the use of a spatial heterodyne method. 
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The presented scheme is the association of a spatial and a temporal heterodyne detection. 

Both temporal and spatial resolution are potentially high. This combination is enabled 
by a wide field measurement performed at one frequency point at a time. The data transfer 
rate bottleneck implies a tradeoff between the number of pixels in the image and the frame 
rate of the detector, which should be guided by the application needs in terms of temporal, 
spectral and spatial resolution. In a few words, combining a spatiotemporal modulation 
and demodulation of a coherent probe light allows one to achieve a sensitive wide field 
detection of a tunable Hertz-resolved spectral component with an array detector. 
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LIST OF FIGURE CAPTIONS 

1. Coherent spectral detection schemes : homodyne (a) and heterodyne (b) optical mix- 
ing. PM : photomixer (square-law detector). SA : spectrum analyzer. Eq object 
field. Ei^o '■ local oscillator field. Si, S2 ■ first and second order object field spectral 
distributions. 

2. Off-axis lensless Fourier configuration for heterodyne holography. L : laser. M : 
mirror. BS : beam splitter. 

3. Measurement of the temporal frequency instrumental response. Camera framerate : 
co's/27r = 8 Hz. Exposure time : = 124 ms. n = 4. Representation of instrumental 
responses for the true image (signal) and the twin image (ghost), in dB. Horizontal 
axis : detuning frequency (culo — ^^L)/(27r), in Hz. Squares : signal (first heterodyne 
term). Circles : ghost (second heterodyne term). 

4. Squared amplitude of the instrumental response defined by eq. [SBl Camera fram- 
erate : uJs/^TT = 8 Hz. Exposure time : Te = 124 ms. n = 4. Representation of 
10 logio[l-^±(^LO — ^^l)!^], in dB. Horizontal axis : detuning frequency (c<jlo— t^L)/(27r), 
in Hz. Dotted line : Continuous line : 

5. Frequency diagram of heterodyne terms spectral components in the case where the 
detuning frequency is set to wlo — ^^l = Acu — us/n. 

6. Images (1024 x 1024 pixels) of the sample for (cjlo — ^^L)/(27r) equal to Hz (a), 400 
Hz (b), 4000 Hz (c) and 8000 Hz (80 ms image exposure time, 4-image demodulation. 
Arbitrary logarithmic scale display. 

7. Traces obtained by summation along columns of FiglHt^a) to (d) intensities. Curves a 
to d correspond to a detuning frequency (culo — ^^L)/(27r) equal to Hz (a), 400 Hz 
(b), 4000 Hz (c) and 8000 Hz (d). Horizontal scale is the image horizontal pixel index. 
Vertical scale is in linear arbitrary units. 

8. Frequency spectra of the light diffused through a suspension of latex particles in brown- 
ian motion. Exposure time is Te = 80 ms. Demodulation is performed with 4 (a), 8 (b). 
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16 (c) and 32 (d) images. Horizontal axis is the detuning frequency (culo — ^^L)/(2vr) 
in kHz. Vertical scale is in linear arbitrary units. The four curves overlap. 

9. Spectra measured with exposure time = 80 ms (a), 20 ms (b), 5 ms (c), and 4-image 
demodulation. Note that curve (a) the same as FigjH] (4 images and 80 ms). Vertical 
scale is in linear arbitrary units. The three curves overlap. 

10. Frequency lineshapes of the light diffused through the cell for different concentrations 
of latex spheres. Exposure time is = 80 ms. 4-image demodulation. Horizontal axis 
is the detuning frequency (c^lo ^ ^l) / (27r) in kHz. Vertical scale is in linear arbitrary 
units. Volumic concentration of latex beads : 2.9 x 10"'^ (a), 1.5 x 10^^ (b), 7.3 x 10~^ 
(c), to 3.6 X 10~^ (d). 
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FIG. 1: Coherent spectral detection schemes : homodyne (a) and heterodyne (b) optical mixing. 
PM : photomixer (square-law detector). SA : spectrum analyzer. Eq object field. -Elo '■ local 
oscillator field, ^i, 5*2 : first and second order object field spectral distributions. 
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FIG. 2: Off-axis lensless Fourier configuration for heterodyne holography. L : laser. M : mirror. 
BS : beam splitter. 
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FIG. 3: Measurement of the temporal frequency instrumental response. Camera framerate : 
u;s/2vr = 8 Hz. Exposure time : Te = 124ms. n = 4. Representation of instrumental responses 
for the true image (signal) and the twin image (ghost), in dB. Horizontal axis : detuning fre- 
quency (wlo — '^L)/(27r), in Hz. Squares : signal (first heterodyne term). Circles : ghost (second 
heterodyne term). 
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FIG. 4: Squared amplitude of the instrumental response defined by eq. [26j Camera framerate : 
iOs/2iT = 8 Hz. Exposure time : Te = 124ms. n = 4. Representation of 101og;^o[l-^±(i^LO — "^l)!^], 
in dB. Horizontal axis : detuning frequency (wlo — '^L)/(27r), in Hz. Dotted line : B^. Continuous 
line : -B_ 
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FIG. 5: Frequency diagram of heterodyne terms spectral components in the case where the detuning 
frequency is set to cjlo — ^j^l = Au; — Los/n. 




FIG. 6: Images (1024 x 1024 pixels) of the sample for (cjlo - t^L)/(27r) equal to Hz (a), 400 
Hz (b), 4000 Hz (c) and 8000 Hz (80 ms image exposure time, 4- image demodulation. Arbitrary 
logarithmic scale display. 
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FIG. 7: Traces obtained by summation along columns of Figl6][^a) to (d) intensities. Curves a to 
d correspond to a detuning frequency (wlo — '^L)/(2vr) equal to Hz (a), 400 Hz (b), 4000 Hz (c) 
and 8000 Hz (d). Horizontal scale is the image horizontal pixel index. Vertical scale is in linear 
arbitrary units. 
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FIG. 8: Frequency spectra of the light diffused through a suspension of latex particles in brownian 
motion. Exposure time is Te = 80 ms. Demodulation is performed with 4 (a), 8 (b), 16 (c) and 32 
(d) images. Horizontal axis is the detuning frequency (wlo — '^L)/(2vr) in kHz. Vertical scale is in 
linear arbitrary units. The four curves overlap. 
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FIG. 9: Spectra measured with exposure time Tg = 80 ms (a), 20 ms (b), 5 ms (c), and 4-image 
demodulation. Note that curve (a) the same as FigJS] (4 images and 80 ms). Vertical scale is in 
linear arbitrary units. The three curves overlap. 




FIG. 10: Frequency lineshapes of the light diffused through the cell for different concentrations 
of latex spheres. Exposure time is Tg = 80 ms. 4-image demodulation. Horizontal axis is the 
detuning frequency (wlo ~ '^L)/(2vr) in kHz. Vertical scale is in linear arbitrary units. Volumic 
concentration of latex beads : 2.9 x lO"^ (a)^ 1.5 x 10"^ (b), 7.3 x lO"'' (c), to 3.6 x 10"^ (d). 
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